Nonequilibrium transport in mesoscopic multi-terminal SNS Josephson junctions 
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We report the results of several nonequilibrium experiments performed on superconduct- 
ing/normal/superconducting (S/N/S) Josephson junctions containing either one or two extra ter- 
minals that connect to normal reservoirs. Currents injected into the junctions from the normal 
reservoirs induce changes in the electron energy distribution function, which can change the prop- 
erties of the junction. A simple experiment performed on a 3-terminal sample demonstrates that 
quasiparticle current and supercurrent can coexist in the normal region of the S/N/S junction. 
When larger voltages are applied to the normal reservoir, the sign of the current-phase relation of 
the junction can be reversed, creating a "7r-junction." We compare quantitatively the maximum 
critical currents obtained in 4-terminal 7r-junctions when the voltages on the normal reservoirs have 
the same or opposite sign with respect to the superconductors. We discuss the challenges involved 
in creating a "Zeeman" 7r-junction with a parallel applied magnetic field and show in detail how 
the orbital effect suppresses the critical current. Finally, when normal current and supercurrent are 
simultaneously present in the junction, the distribution function develops a spatially inhomogeneous 
component that can be interpreted as an effective temperature gradient across the junction, with 
a sign that is controllable by the supercurrent. Taken as a whole, these experiments illustrate the 
richness and complexity of S/N/S Josephson junctions in nonequilibrium situations. 

PACS numbers: 74.50.+r, 73.23.-b, 85. 25. Am, 85.25.Cp 
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I. INTRODUCTION 

When a superconducting metal (S) and a normal metal 
(N) are placed in contact with each other, the proper- 
ties of both metals are modified near the S/N interface. 
This effect, called the superconducting proximity effect, 
was widely studied in the 1960'si* Our microscopic un- 
derstanding of the proximity effect underwent dramatic 
progress in the 1990's as a result of new experiments 
performed on submicron length scales, coupled with the- 
oretical ideas about phase-coherent transport from meso- 
scopic physics. It is now understood that the conven- 
tional proximity effect in S/N systems and the dc Joseph- 
son effect in S/N/S junctions arise from the combina- 
tion of three ingredients: Andreev reflection of electrons 
into holes (and vice versa) at the S/N interface, quantum 
phase coherence of electrons and holes, and time-reversal 
symmetry in the normal metal. Our new understanding 
of the proximity effect in equilibrium situations and in 
linear response transport is demonstrated by a wealth 
of beautiful experiments^ and is summarized in several 
theoretical reviews 

In the past several years, research in S/N systems has 
increasingly focused on nonequilibrium phenomena. Un- 
derstanding nonequilibrium situations is more difficult 
than understanding near-equilibrium situations, because 
the electron energy distribution function in nonequilib- 
rium may be quite different from a Fermi-Dirac func- 
tion. In such situations, the behavior of a specific sample 



may depend critically on the rates of electron-electron or 
electron-phonon scattering. A pioneering work in this 
area was the demonstration by Baselmans et al~ that 
the current-phase relation of a S/N/S Josephson junc- 
tion can be reversed, producing a so-called "7r-junction" . 
This effect is produced by applying a voltage that suit- 
ably modifies the form of the distribution function. 

This paper presents results of several experiments per- 
formed on S/N/S Josephson junctions with extra leads 
connecting the N part of the devices to large normal 
reservoirs. Samples are made from polycrystalline thin 
films of aluminum (S) and silver (N) deposited by thermal 
evaporation. Electrical transport is in the diffusive limit 
- i.e. the electron mean free path is much shorter than all 
other relevant length scales in the problem, including the 
sample length and the phase coherence length. In these 
experiments, the two superconductors are usually at the 
same potential, referred to as ground. Different voltages 
are applied to the normal reservoirs, which in most cases 
cause the distribution function in the structures to devi- 
ate strongly from a Fermi-Dirac distribution. 

Several of the experiments have been analyzed quanti- 
tatively within the framework of the Usadel equations, 4,6 
which are appropriate for S/N samples in the diffusive 
limit. The equilibrium component of the Usadel equa- 
tion is a diffusion equation describing pair correlations 
in N and S. The nonequilibrium, or Keldysh, component 
consists of two coupled Boltzmann equations for the spec- 
tral charge and energy currents. Incorporating inelastic 
scattering into the Keldysh equations involves inserting 
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the appropriate collision integrals; but this procedure has 
so far been followed fully in only a few cases. More- 
over, the effect of inelastic scattering on the equilibrium 
component of the Usadel equation or the proximity ef- 
fect on the collision integrals have never been included 
sclf-consistently to our knowledge. More commonly, re- 
searchers analyzing nonequilibrium phenomena solve ei- 
ther the Keldysh equation without collision integrals, or 
the standard Boltzmann equation with collision integrals 
but without superconducting correlations, depending on 
which aspect of the problem is more important. At the 
end of this paper we compare these approaches as applied 
to the last experiment discussed in the paper. 




FIG. 1: SEM image of sample with two superconducting reser- 
voirs, labeled SI and S2, and normal reservoir labeled N. A 
tunnel probe, labeled TP, consists of thin Al oxidized prior to 
deposition of Ag wire. 

The paper is organized as follows: Section 11 describes 
the sample fabrication and measurement techniques. Sec- 
tion III describes a simple experiment, called the "dan- 
gling arm", involving a 3-terminal S/N/S device with a 
single extra lead to a normal reservoir. The dangling arm 
experiment was first reported by Shaikhaidarov et al.. 7 
We include it here because it provides a clear demon- 
stration of the superposition of quasiparticle current and 
supercurrent in a S/N/S junction, an essential result for 
the remainder of the paper. Section IV describes the w- 
j unction experiment in 3- and 4-terminal devices. The 
4-terminal sample allows a direct comparison of the situ- 
ations present in the 3-terminal 7r-junction£ and the orig- 
inal 4-terminal 7r-junction of Basclmans et aZ.— Section V 
discusses the behavior of the critical supercurrent as a 
function of magnetic field applied parallel to the plane 
of the sample, and shows the difficulty involved in try- 
ing to achieve a 7r-junction by Zeeman splitting of the 
conduction electrons .^£1 The theoretical calculation rel- 
evant to this geometry is given in the appendix. Section 
VI discusses an experiment in which supercurrent and 
quasiparticle current are independently controlled in a 
3-terminal S/N/S junction, leading to an effective tem- 
perature gradient across the junction^ The local distri- 



bution function is measured by a tunnel probe near one 
of the S/N interfaces. The discussion provides informa- 
tion that was not included in our previous report on this 
experiment. 1 ^ Together these experiments demonstrate 
the richness of phenomena present in S/N/S Josephson 
junctions under nonequilibrium conditions. 



II. EXPERIMENTAL TECHNIQUES 
A. Fabrication 

All samples in this work were fabricated using e-beam 
lithography. A bilayer of resist was deposited onto an un- 
doped Si wafer covered only with its native oxide layer. 
The bilayer was formed by first depositing a copolymer 
P (MM A/MA A), followed by a second layer of PMMA. 
The bilayer was exposed by 35-keV electrons and then 
developed to make a mask for evaporation. With the 
resist bilayer, it is possible to fabricate undercuts in 
the mask, allowing angled evaporation techniques to be 
used. 13 Therefore, multiple layers of different metals (ei- 
ther 99.99% purity Al or 99.9999% purity Ag) were se- 
quentially deposited without breaking vacuum. 

These techniques were used to prepare the sample 
shown in Fig. [TJ To create the tunnel probe (TP), 30 
nm of Al was deposited while the sample was tilted 45 
degrees, creating an actual thickness of about 21 nm of 
Al on the surface. Next, a mixture of 90%Ar-10% O2 
gas was leaked into the vacuum chamber to a pressure 
of 60 Torr. After 4 minutes, the chamber was evacuated 
again, in preparation for the following depositions: For 
the silver wires; labeled R±, R2, and Rn', 30 nm of Ag 
was deposited with the plane of the wafer perpendicu- 
lar to the evaporation source. For the superconducting 
reservoirs, SI and S2, the sample was tilted 45 degrees 
and rotated 180 degrees in order to deposit 90 nm of 
Al (for a 60 nm thickness). Finally, the sample was ro- 
tated another 140 degrees in preparation for a final, thick 
layer of Ag to be deposited over the normal reservoir, N. 
The sample in Fig. [S] followed a similar procedure, except 
foregoing the first Al deposition and oxidation steps. 



B. Experimental setup 

Samples were measured inside the mixing chamber of 
a top loading dilution refrigerator. All electrical leads 
to the sample passed through commercial LC 7r-filters at 
the top of the cryostat and cold RC filters in the cryostat 
consisting of 2.2 kSl resistors in series and 1 nF capacitors 
coupled to ground. 

Current-voltage characteristics (I-V curves) were ob- 
tained through 4-probe measurements across the sample. 
The current was swept using a triangle wave and several 
cycles were collected and averaged together. Measure- 
ments of dl/dV were obtained by adding a slow (~ 1 
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mHz) triangle wave pattern to the sine output of a lock- 
in amplifier. The lock-in amplifier was operated at low 
frequencies (less than 100 Hz) to allow for extrapolation 
of the system response to zero frequency. Both the in- 
phase and out-of-phase components of the signal were 
recorded and utilized in the analysis. 



III. DANGLING ARM EXPERIMENT 

The dangling arm experiment was first proposed in 
the Ph.D. thesis of S. Gueronjii although a related ge- 
ometry was discussed by Volkov two years earlier^ The 
experiment is performed on a 3-terminal S/N/S Joseph- 
son junction sample similar to the one shown in Fig.[TJ in 
which the tunnel probe in the lower left was unused. We 
label the three terminals of the sample Si, 82, and N, and 
the resistances of the three arms Ri, R2, and Rn- (We 
neglect for the moment the variation of these resistances 
due to proximity effect.) One measures the resistance 
from N to Si, while leaving S2 open (dangling). Naively, 
one might expect the measured resistance between N and 
Si to be equal to Rn + Ri- That result would imply that 
the current travels directly from N to Si, which in turn 
implies that Si and S2 are at different voltages. Given 
that Si and S2 are coupled by the Josephson effect, the 
relative phase (f> between Si and S2 then accumulates at 
the Josephson frequency, d(j)/dt = leVnjh. If, however, 
the injected current / splits into a piece 1\ through R\ 
and a piece I2 through R2, such that I\R\ = /2-R2, then 
Si and S2 will be at the same potential. To avoid having 
a net current flowing into the dangling arm, the sample 
must then provide a supercurrent 1$ from S2 to Si that 
exactly cancels the quasiparticle current 12- In this sce- 
nario, R\ and R2 are effectively acting in parallel, and the 
measured resistance will be Rp = Rn + R1R2/ (Ri + R2)- 

Figure [2] shows the 2-terminal I — V curve taken at 
T = 51 mK from a sample similar to the one shown 
in Fig. [TJ with nominal resistance values R\ = 7.0 ft, 
R 2 = 7.0 fl, and R N = 16.9 fl. The inset shows dV/dl 
vs. /, providing a clearer view of the effective resis- 
tance. Either plot shows that the resistance is about 
20.7 f2 when the applied current is less than about 0.94 
/iA. This resistance is very close to the nominal value of 
Rp = 20.4 fi. When the current exceeds 0.94 /LtA, the re- 
sistance increases to the value 24.6 f2, which is very close 
to Rn + Ri — 23.9 il. (Resistance differences less than an 
Ohm are attributed to the finite size of the "T-junction" 
in the middle of the sample.) These data confirm the 
idea outlined in the previous paragraph, that supercur- 
rent and quasiparticle current can coexist in the normal 
region of a S/N/S Josephson junction. 

The transition at 1^ where the resistance increases 
to Rn + R± occurs when the supercurrent across the 
S/N/S Josephson junction exceeds the S/N/S critical 
current, I^ NS . However, since only the fraction I2 — 
IRi/(Ri + R2) ~ 1/2 of the injected current must be 
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FIG. 2: Voltage versus current measured between reservoir N 
and Si with S2 floating, at T — 51 mK. Dotted lines represent 
slopes of 20.7 and 24.6 fi, which correspond to the resistances 
Rp and Rn + Ri, respectively. Inset: Differential resistance 
vs. current under similar conditions, showing agreement be- 
tween the two measurement techniques. 



cancelled by the supercurrent, one should expect that 
= I^ NS (R 1 + r 2 )/r 1 . The data in Figure [3] show 
that this expectation is fulfilled at relatively high tem- 
peratures, but that at lower temperature falls well 
below this value. 

Two reasons for the small values of at low tem- 
perature were given by Shaikhaidarov et at. 7 Those 
authors solved the Usadel equation analytically in the 
limit where the S/N interfaces have high resistance (poor 
transparency), so that proximity effects are small and 
the Usadel equation can be linearized. They pointed out 
that Ic is suppressed below its equilibrium value due 
to the applied voltage at N, a result we will reinforce 
below. They also argued that the phase-dependence of 
the resistances R%, i?2, and Rn due to proximity effect 
causes the measured value of I^ s to be smaller than the 
nominal value I? s = I^ NS (R 1 + R 2 )/Ri- 

We believe that the effect related to the phase depen- 
dence of the resistances is small and the relative decrease 
in I^ s at low temperature is due predominantly to the 
decrease in I^ NS as a function of U. This effect is demon- 
strated graphically in the inset to Fig. [31 There the crit- 
ical current of the Josephson junction, I^ NS : multiplied 
by the constant ratio (i?i + R2)/ R\, is plotted as a func- 
tion of the voltage U applied to the normal reservoir. 
As can be seen in the inset, if NS decreases rapidly as 
a function of U. The straight line through the origin in 
the inset represents the current injected into the sam- 
ple from the N reservoir, U/Rp, where the resistances 
are evaluated at phase difference 7r/2 between Si and S 2 . 
The intersection of the two curves shows the value of the 
dangling arm critical current (ordinate) at the criti- 
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FIG. 3: Critical current measured between N and Si (•), 
between N and S2 (a), and between Si and S2 (□) - the 
latter multiplied by the ratio (Ri + -R2V-R1 ~ 2 - for different 
temperatures. The three data sets are in close agreement at 
temperatures above about 250 mK. Inset: Graphical approach 
to calculation of low-temperature critical current between N 
and Si. The dots are measurements of the critical current 
across S1-S2, again multiplied by (Ri + R2) / Ri, as a function 
of applied voltage U between N and Si. The critical current 
decreases rapidly with increasing U. The line through the 
origin represents the injected current from N. The intersection 
gives the critical current 1^ at the critical value of U. Note 
that all critical current values in the inset are 15 — 20% larger 
than in the main panel, due to a small magnetic field B = 125 
G present when the latter data were obtained. 

cal voltage U^ s (abscissa). The figure demonstrates the 
large reduction in S/N/S critical current due to the ap- 
plied voltage U, which explains why is much smaller 
than I^ NS (Ri + R 2 )/R\ at low temperature. At high 
temperatures T > eU/ks, the relative reduction is less 
significant due to two reasons: First, increasing the tem- 
perature decreases the critical current if NS , and thereby 
also U^ s . Moreover, to observe a sizable reduction in 
l! NS {U), \eU\ has to exceed k B T. 

IV. S/N/S NONEQUILIBRIUM tt-JUNCTION 

A. Three-terminal 7r-junction 

Figure [3] shows, not surprisingly, that the critical cur- 
rent I c of an S/N/S Josephson junction decreases when 
quasiparticle current is injected into the junction from 
a normal reservoir. Indeed, if the only effect of the in- 
jected current were to heat the electrons in the junc- 
tion, then one would expect the critical current to con- 
tinue decreasing monotonically as a function of the ap- 
plied voltage £/J£ That this is not the case represents 



a major discovery in nonequilibrium superconductivity 
by Baselmans et al& in 2000. Those authors showed 
that I c first decreases as a function of U, but then in- 
creases again at higher U. The explanatio n 1 ^ 18 for this 
counter-intuitive result consists of two pieces. First, one 
can view the supercurrent in the sample as arising from 
the continuous spectrum of Andreev bound states in the 
normal metal ) 19 i 20 which carry supercurrent in either di- 
rection, depending on their energy. Second, in the pres- 
ence of the applied voltage U the electron distribution 
function in the junction is not a hot Fermi-Dirac distri- 
bution, but is closer to a two-step distribution - as long as 
the short sample length does not allow electron thermal- 
ization within the sampled The two-step distribution 
function preferentially populates the minority of Andreev 
bound states that carry supercurrent in the direction op- 
posite to the majority, hence it reverses the current-phase 
relation in the junction i 17 i 18 Such a Josephson junction 
is called a "7r-junction" , because the energy-phase and 
current-phase relations arc shifted by tt relative to those 
of standard Josephson junctions. 



The original 7r-junction experiment of Baselmans et al. 
was performed in a 4-terminal sample, where voltages of 
opposite sign were applied to the two normal reservoirs. 
Later, Huang et al. 8 demonstrated that a 7r-junction can 
also be obtained in a 3-terminal geometry with a single 
normal reservoir, a result predicted by van Wees et al^ 
10 years earlier. 



Figure S^, shows results of a 3-terminal 7r-junction ex- 
periment performed on a sample similar to the one in 
Fig. [TJ where the tunnel probe in the lower left portion 
of the figure is not used. We measure the I — V curve of 
the S/N/S Josephson junction using a 4-probe current- 
bias measurement, while a dc voltage is simultaneously 
applied to the normal reservoir via a battery-powered 
floating circuit. Figure S^l shows a series of I — V curves 
at different values of the voltage U applied to the nor- 
mal reservoir. Figure [Dd shows the critical current I c vs. 
U . Notice that I c initially decreases with an increasing 
U, as shown in the inset to Fig. [3J But as U increases 
further, I c reaches a minimum value (indistinguishable 
from zero in this experiment^) at U = U c w 34 fiV, then 
grows again to reach a second maximum at U ~ 63 /iV. 
The minimum in I c separates the standard Josephson 
junction behavior at low values of U from the 7r-junction 
behavior at higher U. If instead of plotting the critical 
current I c (which is by definition a positive quantity) one 
were to plot the supercurrent I s at a fixed phase differ- 
ence (j) = 7r/2 across the junction then the graph would 
show a smooth curve passing through zero at U = U c , 
reaching a local minimum at U ~ 63 /uV and gradually 
returning to zero at large U. 
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FIG. 4: a) Voltage vs. current across the S/N/S junction for 
selected voltages, U , applied to the normal reservoir. Graphs 
for different U are offset for clarity, with (7 = 0, 17, 29, 35, 
41, 63, 92, and 114 //V from bottom to top. The hysteresis in 
the (7 = data is probably due to heating of the Ag wire in 
the normal state, b) Critical current vs. U. 



B. Comparison of four-terminal 7r-junctions with 
symmetric and antisymmetric bias 

The physical explanation of the 7r-junction in the 3- 
tcrminal sample is the same as in the 4-terminal sample, 
with the differences arising only from the distribution 
functions. Figure O shows a schematic drawing of the 
distribution function f(E) along a path from a reser- 
voir N to S for both 4-terminal and 3-terminal samples 
for U » ksT, assuming weak electron-electron inter- 
actions in the N wire, and neglecting the proximity cor- 
rections. Notice in figures c and d that f(E) consists 
of a double-step function, the step height within the en- 
ergy range — elf to elf changing with the location along 
the wire. As we will show in the next section, the even 
(in energy) part of f(E) has no effect on the magnitude 
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FIG. 5: a) Depiction of electron flow in four-terminal configu- 
ration, b) Depiction of electron flow in three-terminal config- 
uration, c) Schematic representation of the distribution func- 
tion on the path between a normal (N) and superconducting 
(S) terminal in the structure (a) under high bias U » ksT. 
d) Distribution function in the T-structure (b) under similar 
conditions. Due to Andreev reflection f(e) is discontinuous 
at the N-S interface as explained in the text. 



of the supercurrent (in the absence of electron-electron 
interactions) , suggesting that the voltage-dependent crit- 
ical current, I C (U), would be identical in 3-terminal and 
4-terminal samples with identical dimensions and resis- 
tances. However, there are three reasons why this is not 
quite true: First, Joule heating is more prevalent in the 
4-terminal device, which rounds the distribution func- 
tions more than in the 3-terminal device. Second, the 
spectral supercurrent density, jE(E), evaluated at the 
junction point will be slightly smaller in the 4-terminal 
sample than in the 3-terminal sample due to the pres- 
ence of the additional arm connecting the sample to a 
normal reservoir^ Finally, f(E) will be slightly more 
rounded in the 4-terminal sample due to the increased 
phase space available for electron-electron interactions. 
Roughly speaking, the rate of e-e interactions at a given 
energy is proportional to f(E)[\ — f(E)], which is max- 
imized when f(E) = 1/2. Each of these effects serve 
to increase I c in the 3-terminal geometry relative to the 
4-terminal geometry. 

It is not practical to compare critical currents from 
two different samples, since they will never have identi- 
cal dimensions nor electrical resistances. Instead, it was 
proposed in Ref. 8 to compare the critical currents in a 
single 4-terminal sample under conditions of symmetric 
and antisymmetric voltage bias of the two normal reser- 
voirs. Figure [5] shows the sample we fabricated for this 
experiment, with the superconductive reservoirs labeled 
SI and S2 and the normal reservoirs labeled Nl and N2. 
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By applying a positive potential U to Nl and a nega- 
tive potential — U to N2, one reproduces the experiment 
performed by Baselmans et al. We call this situation 
antisymmetric bias, since the two applied voltages differ 
by a negative sign. In this case the quasiparticle cur- 
rent overlaps with the supercurrent only at the crossing 
point of the sample where the electrostatic potential is 
equal to zero. In contrast, applying the identical voltage 
U on both Nl and N2 (with ground defined at one of 
the superconducting electrodes), called symmetric bias, 
will produce a situation mimicking that in the 3-terminal 
experiment of Huang et al& Notice that by mimicking a 
3-terminal sample with a 4-terminal sample, geometrical 
differences between the two experiments are eliminated, 
so any observed difference in the critical currents will be 
due either to e-e interactions or to Joule heating. 
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FIG. 6: SEM image of S/N/S Josephson junction. The *x'- 
shape is deposited Ag that connects to (difficult to see) Al 
reservoirs above and Ag reservoirs below patterned by angled 
evaporation. The feature in the Ag wire near N2 is likely due 
to a near burn in the sample. 

The preceding description of symmetric and antisym- 
metric biases holds strictly only if the resistances of the 
two lower arms are identical. Otherwise, application of 
antisymmetric bias will result in a nonzero potential at 
the cross and some quasiparticle current will flow into the 
superconducting reservoirs. In that case, f(E) will take 
a form intermediate between those depicted in Figures 
[SJ; and d, which decreases the measurement contrast be- 
tween the two biases. In our experiments we took care 
to measure the resistances of all the arms and to ensure 
that the voltages at the two normal reservoirs were indeed 
equal (for symmetric bias) or opposite (for antisymmetric 
bias). 

Figures [Th. and [7b show I —V curves measured across 
the S/N/S junction at T = 170 mK, for several different 
values of U. Curves with increasing values of U are offset 
upward for clarity. Figure [7^, shows the data for anti- 
symmetric bias while Fig. [7)d shows symmetric bias. The 
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FIG. 7: Data showing voltage drops across segments of wire 
either in antisymmetric or symmetric arrangements while cur- 
rent flows from SI to S2. Each line (offset for clarity) corre- 
sponds to a different value of U. a) Voltage across SI to S2 
for the antisymmetric measurement. Applied voltages U are 
from the bottom: 19, 28, 38, 52, and 71 fiV b) Voltage across 

51 to S2 for the symmetric measurement. Applied U: 17, 25, 
37, 49, 72, and 131 /iV. c) Resistance across Nl to N2 for 
the antisymmetric measurements, d) Resistance across N2 to 

52 for the symmetric measurements, taken from voltage mea- 
surements in which a constant resistance was subtracted from 
the graph. 



data follow the same trend observed in Fig. [4] namely, the 
critical current first decreases with increasing U, then in- 
creases again before finally disappearing altogether. Fig- 
ure^ shows this critical supercurrent as a function of U. 
In this figure we have plotted the critical current as neg- 
ative in the range of voltages after I c disappears initially, 
to signify that the junction is in the 7r-state as discussed 
earlier. 

The transition from the 0-state to the 7r-state can be 
confirmed directly by experiment, 5 without recourse to 
the theoretical explanation. The resistances of the Ag 
arms of the sample vary with the phase <f> due to the 
proximity effect^ The phase cf>, in turn, varies between 
±7r/2 as a function of the supercurrent Is passing be- 
tween Si and S2; hence, one observes a variation of the 
resistances as a function of Is- This effect is shown in 
Figs. \7]p and d, in which the resistances between Ni and 
N2 (Ni and S2) were measured versus Is for the antisym- 
metric and symmetric bias configurations, respectively. 
Each curve in the lower two figures has the same value 
of U as the corresponding I — V curve in the upper fig- 
ures. One can see that proximity effect induces a local 
minimum in the resistance at I s = when the junction is 
in the 0-state, because that is where = 0. In contrast, 
the resistance exhibits a local maximum in the resistance 
at / s = when the junction is in the 7r-state, because 
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FIG. 8: (color online) Critical current of a 4-terminal S/N/S 
Josephson junction versus voltage U applied to the normal 
reservoirs. The voltages are applied either antisymmetrically 
(■) or symmetrically (a) to the two reservoirs. Solid lines 
represent simultaneous best fits to data at different tempera- 
tures. Fitting methods are discussed in the text for data taken 
at bath temperatures (a) 35 mK and (b) 170 mK. The dashed 
line represents the best fit when Joule heating is excluded. 



4> = ir. Interestingly, the top curve in Fig. [TJi shows that 
at large enough values of U, the system returns to the 
0-state since the resistance again shows a local minimum 
at I s = 0. This high-[/ transition from the 7r-state back 
to the 0-state was not visible in the S/N/S I —V curves. 

Figure [5] shows the behavior I c vs. U at two different 
temperatures. The squares represent antisymmetric bias 
while triangles represent symmetric bias. Both bias con- 
figurations appear similar in that the samples cross to 
the 7r-state at nearly the same value of U. It should be 
noted, however, that the maximum 7r-current is larger for 
symmetric bias than for antisymmetric bias. That result 
is consistent with the qualitative arguments made above. 
In the next section we present a quantitative analysis of 
the results. 



C. Calculation of the Critical Current in an S/N/S 
Josephson Junction 

The amount of supercurrent passing through an S/N/S 
Josephson junction may be calculated byH 

<7 N A 



Is 



dE[l-2f(E)]j E (E) (la) 



= <j n A / dE f L (E)j E (E) 



(lb) 



where ctjv and A are the conductivity and cross-sectional 
area of the normal metal, respectively. f(E) is the dis- 
tribution function within the normal wire, and /l(E) = 
f(—E) — ,f(E) is the antisymmetric component of f(E) 
with respect to the potential of the superconductors. The 
spectral supercurrent density, je(E), is an odd function 
of energy, and describes the amount of supercurrent at 
a given energy travelling between superconductors with 
relative phase difference <j>. In the samples considered 
in the present work, it is generally sufficient to calculate 
the supercurrent using the the distribution function at 
the crossing point of the wires. 

To determine ]e{E), we solve the Usadel equation nu- 
merically using the known physical dimensions and elec- 
trical resistances of the various wire segments of the sam- 
ple. We then look for consistency with the measured 
temperature dependence of the equilibrium critical cur- 
rent, I C (T) shown in Fig. [9] The je{E) used to fit these 
data, evaluated at <j> = n/2, is shown in the inset of Fig. [9] 
Since the length L of the junction is much longer than the 
superconducting coherence length £ s of the S electrodes 
for all the samples studied in this work, the damped os- 
cillations in ]e{E) occur on an energy scale given by the 
Thouless energy, Exh = frD/L 2 , where D is the diffusion 
constant in the wire. Exh characterizes the temperature 
scale over which the equilibrium critical current drops to 
zero, and also determines the voltage scale U needed to 
create a nonequilibrium 7r-j unction. The transition from 
the 0-state to the 7r-state occurs at ell w 8Exh- The fit 
shown in Fig. [9] was obtained with Etk =4.11 /ieV. 

Next we calculate f(E) in the nonequilibrium situation 
with antisymmetric bias, i.e. with voltages U and —U ap- 
plied to reservoirs Ni and N2. Because we are interested 
in the distribution function far from the superconducting 
reservoirs, we consider f(E) using the Boltzmann equa- 
tion. Let us first ignore the supercurrent and proximity 
effects, although inclusion of those effects will be dis- 
cussed in detail in Sec. IVII In a reservoir at voltage U, 
f(E) is a Fermi-Dirac distribution displaced by energy 
eU, f(E) = f FD (E-eU) = [exp((£ - eU)/k B T) + 1]" 1 . 
In the experiment with antisymmetric bias, we then have 
f(E) = f FD {E - eU) at reservoir Ni and f(E) = 
Jfd(E + eU) at reservoir N 2 . If we neglect inelastic 
electron scattering, then in the middle of the wire (at 
the intersection of the cross) f(E) has the double-step 
shape: 

f(E) = l[f FD (E + eU) + f FD (E-eU)} (2) 
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FIG. 9: Critical current I c at several temperatures for the 
sample shown in Fig. [6] The line is the best fit by solving 
equation llbl with a Fermi-Dirac distribution function. Inset: 
Solution for spectral supercurrent used to fit data. 



In fact, the odd part of f(E) is the same everywhere in 
the Ag wire: 



h(E) 



tanh 



E-eU 
2k B T 



tanh 



E + eU 
2k B T 



(3) 



This conclusion holds also in the presence of the proxim- 
ity effect. At energies \E\ < ell, the even part, fr(E), 
varies linearly with distance between the two reservoirs 
(the lower two arms of the sample), but is zero every- 
where along the direct path connecting the two super- 
conductors (in the ideal case where the resistances of the 
two lower arms are equal). 

To calculate f(E) in the experiment with symmetric 
bias, we need the boundary conditions at the interfaces 
between the normal wires and the superconducting reser- 
voirs. For energies below the superconducting gap, A, 
these conditions are fr = and O/l/Ox = 0, where 
f T (E) = 1 - f(E) - f{-E). These boundary condi- 
tions assume high-transparency interfaces, no charge im- 
balance in the superconductors, and no heat transport 
into the superconductors^ (Note that Jl(E) is discon- 
tinuous at the N-S interface for energies below the gap, 
and returns to the standard form tanh(E /2ksT) in the 
S electrodes.) The solution for f(E) at the N/S inter- 
face is identical to Eq. ||2J), but the symmetric compo- 
nent /t(E) is nonzero elsewhere in the wire. Notice that 
the odd component of the distribution function, f^(E), 
is identical in the two cases everywhere in the sample. 
The proximity effect induces a small feature in /l(E) 
discussed in Sec. IVI1 but it is zero at the crossing point 
in the middle of the sample. 

Calculation of f(E) in the realistic situation requires 
consideration of electron-electron interactions in the Ag 



wire. (The electron-phonon interaction, in contrast, is 
much weaker, and need be considered only in the mas- 
sive normal reservoirs. See the discussion below.) To 
incorporate electron-electron interactions, we solved the 
Boltzmann equation in the wire numerically, following 
previous work by Pierrej 26 ' 27 The results of this numeri- 
cal calculation of f(E) in the situations with either sym- 
metric or antisymmetric bias were extremely similar. In- 
deed, the slight additional rounding of f(E) in the an- 
tisymmetric case could not account for the differences 
observed in the experiment, shown in Fig.[5J 

To account for the difference in the observed I C (U) 
between the two experiments, we next considered the ef- 
fect of Joule heating on the temperatures of the normal 
reservoirs. (Due to Andreev reflection at the N/S inter- 
faces, there is no heat transport into the superconduct- 
ing reservoirs.) Although we intentionally fabricated the 
normal reservoirs much thicker than the wires, this was 
not enough to eliminate the effects of Joule heating al- 
together. The heat current in a reservoir at a distance r 
from the juncture with the wire is given by: 



3 



Q 



-£<jT\7T: 



P 
Ifrt 



(4) 



where P = I 2 R is the total power dissipated in the wire, 
a and t — 310 nm are the conductivity and the thickness 
of the reservoir, respectively, and £ = 7r 2 /3(fcs/e) 2 = 
2.44 x 10~ 8 V 2 /K 2 is the Lorenz number. (We neglect 
the small additional Joule heat generated in the reservoirs 
themselves.) The spreading angle 6 « n if we consider 
the combination of the two normal reservoirs shown in 
Fig. El 

Using ([4| as a boundary condition, one can find an 
effective temperature at the wire-reservoir interface equal 
tc 



,28 



Toff = VT 2 + b 2 U 2 



(5) 



The temperature far away in the normal reservoir is as- 
sumed to be T, the bath temperature. The factor b is 
given by22 



2 = R a n 
6»£i? n r ' 



(6) 



where i?n = 1/ (at) is the sheet resistance of the normal 
reservoir, r$ w the wire width, and r\ is the distance over 
which the electrons in the reservoir thermalize to the bath 
temperature via electron-phonon scattering. The param- 
eter b varies inversely with the thickness of the metal 
reservoir and the electrical resistance encountered by the 
quasiparticle current in the wire. Because the voltage 
drop U in the antisymmetric bias situation occurs en- 
tirely between a normal reservoir and the crossing point, 
the resistance R is smaller than in the symmetric bias 
situation where U drops fully from the N reservoirs to 
the N/S interfaces. The larger current in the former case 
causes more Joule heating, and hence a larger reservoir 
temperature. For our sample, the values of b needed to 
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fit the data (see solid lines in Figs. [8^ and [8)3) are 2.7 
and 3.2 K/mV, respectively, for the symmetric and anti- 
symmetric bias experiments. Their ratio of 1.2 matches 
the ratio calculated from the sample parameters. Their 
magnitudes, however, are nearly three times larger than 
what we calculate based on the total reservoir thickness. 
The experiment seems to suggest that heat was trapped 
in the 35-nm Al layer at the bottom of the reservoirs, 
rather than immediately spreading throughout the whole 
reservoir thickness^ 



ferromagnet. In the standard elemental ferromagnets, 
E ex is large (ss 0.1 meV), hence £f is extremely short - 
on the order of 1 nm. Control of sample thickness unifor- 
mity at this scale is difficult, hence several groups have 
used dilute ferromagnetic alloys, with reduced values of 
E ex , to increase £p. The advantage of the "Zeeman" 
7r-junction is that it is fully tunable by the field. The 
disadvantage is that the sample must be thin enough to 
minimize the effects of orbital pair-breaking in both the 
superconducting electrodes and in the normal part of the 
junction. 



V. APPLICATION OF A PARALLEL 
MAGNETIC FIELD, AND THE "ZEEMAN" 
7r- JUNCTION 

There is a long history of applying magnetic fields 
perpendicular to the direction of current flow in super- 
conductor/insulator/superconductor (S/I/S) Josephson 
junctions, to observe the famous Fraunhofer pattern in 
the critical current. In S/N/S junctions, the Fraunhofer 
pattern is observed only in wide junctions, whereas nar- 
row junctions exhibit a monotonic decrease of the crit- 
ical current with field due to the orbital pair-breaking 
effect.^ 

In this section we discuss the effect of a magnetic field 
applied parallel, rather than perpendicular, to the cur- 
rent direction. In this geometry there should never be a 
Fraunhofer pattern. And because the samples are thin 
films, one expects the orbital pair-breaking effect to be 
much weaker than for a field applied perpendicular to the 
plane. In the case of an extremely thin sample the Zee- 
man (spin) effect should dominate over the orbital effect 
of the field. 

The effect of Zeeman splitting on an S/N/S Joseph- 
son junction was studied theoretically in 2000 by Yip 10 
and by Heikkila et al.S- Their idea is that the electronic 
structure of a normal metal in a large applied magnetic 
field resembles that of a weak ferromagnet, in that the 
up and down spin bands are displaced by the Zeeman 
energy. They also showed how the Zeeman-split junc- 
tion behaves analogously to the nonequilibrium S/N/S 
junction, with the Zeeman energy playing the role of the 
voltage in Eqs. (|lb[) and Josephson junctions made 
with real ferromagnetic materials (S/F/S junctions) are 
the subject of intense current interest, as they can also 
show 7r-junction behavior^ Unlike the 7r-junctions dis- 
cussed earlier in this paper, however, the 7r-junctions in 
S/F/S systems occur in equilibrium. They appear only in 
particular ranges of the F-layer thickness, due to spatial 
oscillations in the superconducting pair correlations in- 
duced in the F metal near the F-S interface by proximity 
effect. Those oscillations, in turn, arise from the different 
Fermi wavevectors of the spin- up and spin-down electrons 
in the F metal. In diffusive S/F/S junctions, the sign of 
the coupling between the two superconductors oscillates 
over a distance scale £f = (fiD / 'E ex ) x l 2 , where D is the 
diffusion constant and E ex is the exchange energy in the 
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FIG. 10: Critical current across S/N/S junction as a func- 
tion of external magnetic field applied parallel to the cur- 
rent direction. In the absence of orbital pair-breaking effects, 
the transition into the 7r-state would be expected near 0.35 
T. Markers indicate experimental results and solid line the 
scaling (J7J obtained from the Usadel equation. Inset: Criti- 
cal current versus applied voltage for same sample, showing 
transition to the 7r-state at U — 20 fj,V. 

Figure [TU] shows a plot of I c vs. B in an S/N/S sam- 
ple whose normal part had length L = 1.4 /im, width 
w = 50 nm and thickness t = 33 nm. The criti- 
cal current decreases monotonically to zero, over a field 
scale of 0.1 T. This result might appear surprising at 
first glance: At a field B = 0.1 T, the magnetic flux en- 
closed in the cross-section of the wire perpendicular to 
the field is only $ w 1.6 • lO" 16 Tin 2 = 0.08$ , where 
$0 = h/2e is the superconducting flux quantum. Fur- 
thermore, separate tests of the Al banks confirm that 
they remain superconducting to fields of order 0.85 T. 

A quantitative understanding of the data in Figure fTOl 
can be obtained from a solution to the Usadel equation. 
The analysis discussed in Appendix [X] shows how a par- 
allel magnetic field can be absorbed into a spin-flip rate 
r s f in the equations. This allows us to apply the scaling 
I C (B)/I C (B = 0) w cxp(-0.145r sf /£; T/l ) for the zero- 
temperature supercurrent of an S/N/S junction found in 
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Ref. [33j and find 

I C (B)/I C (B = 0) « e 



-(S/Bi) 2 



B x w 6.43 



f 2 



eLwt 



0.10 T. 



(7) 
(8) 



Our numerical calculations confirm that this scaling also 
applies in our multi-probe experimental geometry. This 
prediction is in a good agreement with the experiment, 
as seen in Fig. [TDJ 

In the limit w t, the characteristic field scale B\ 
varies as &o/Lt, rather than the more intuitive result 
$o / wt we might expect based on the cross-sectional area 
of the wire perpendicular to the field. The physical ex- 
planation for this result was given by Scheer et ali^ in 
a paper discussing universal conductance fluctuations as 
a function of parallel field in normal metal wires. As an 
electron travels down the length of a long diffusive wire, 
its trajectory circles the cross-section of the wire many 
times - on order N m (L/w) 2 . Because diffusive mo- 
tion can be either clockwise or counterclockwise as seen 
looking down the wire, the standard deviation in net flux 
and accumulated phase between different trajectories is 
approximately proportional to Bwty/W = BLt, which 
gives the scaling for dephasing. 

It is instructive to ask what constraints on the sam- 
ple geometry would have to be met to enable observa- 
tion of the Zeeman 7r-junction. We estimate the Thou- 
less energy of the sample discussed in Fig. [10] to be 
Etii ~ 2.5 /xeV both from the temperature dependence 
of I c (not shown), and from the voltage- induced tran- 
sition to the 7r-state at U c — 20 /zV (inset to Fig. [TU)) . 
According to theory, the Zeeman tt-j unction should oc- 
cur when gfi B B « 16E Th ^2., or B = 0.35 T. Attempts 
to make thinner samples in order to increase the field 
scale Bi in Eq. Q were unsuccessful, due to the ten- 
dency of very thin Ag films to ball up. According to 
the theory, much thinner films, with t/L of the order of 
0.2gnB/eD ~ 0.001, will be required to enable observa- 
tion of the Zeeman 7r-junction. 



VI. ENGINEERING THE DISTRIBUTION 
FUNCTION 

The discussion in Section IV. B. implied that the 3- 
tcrminal and 4-tcrminal 7r-junctions are similar, with 
only minor differences due to a slight decrease in the 
phase space available for electron-electron interactions in 
the 3-terminal case. But that oversimplified discussion 
misses some important physics. Heikkila et alJ^ showed 
that the superposition of quasiparticle current and super- 
current in the horizontal wire in the 3-terminal sample 
induces a change 5f(E) in the distribution function at 
energies of order Eth- The new feature is antisymmetric 
in space and energy (see Fig.[l6]for the theoretical predic- 
tion and our experimental results which follow each other 
nicely) , and can be interpreted as a gradient in the effec- 
tive electron temperature across the S/N/S junction. For 



this reason, the result was dubbed a "Peltier-like effect." 
Although a tiny cooling effect does occur, observing it 
in a real electron temperature would require a slightly 
modified experimental setup^ In the present case, one 
should view this effect mostly as a redistribution of the 
Joule heat generated in the sample by the applied bias 
U. 

In Sec. lIV Cl it was discussed how the distribution func- 
tions behave in the absence of proximity effects and su- 
percurrent. Including these effects, but ignoring inelastic 
scattering, results in the kinetic equations^ 



— = 0, JT = D T ( X )—+j E f L + T(x)—: 



djL „ . _ n , ,df L . df T 
— = 0, 3L = D L { X )—+ jE f T -T{x) — ; 



(9a) 



(9b) 



where ]t{E) and Jl{E) are the spectral charge and en- 
ergy currents, respectively. The energy-dependent coef- 
ficients Dt, Dl, je, and T can be calculated from the 
Usadel equatior&2&, and vary with the superconducting 
phase difference <fi between Si and S2. In general, these 
equations must be solved numerically; however, they can 
be solved analytically by ignoring the energy dependence 
in Dt and Dl and neglecting the T terms. One can show 
that, in the presence of both the applied voltage U and 
a nonzero supercurrent between Si and S2, /l along the 
horizontal wire connecting the two superconductors con- 
tains a spatially antisymmetric contribution proportional 
to Je{E). In the exact numerical solution to Eqs. ©, the 
feature is distorted due to the rapid evolution of the dif- 
fusion coefficients Dt and Dl near the N/S interfaces. 11 

The antisymmetric feature in f(E) can be measured by 
performing tunneling spectroscopy with a local supercon- 
ducting tunnel probe, which has been demonstrated by 
Pothier et al. 21 to reveal detailed information about f(E) 
in a metal under nonequilibrium conditions. In our case, 
the local probe must be placed close to the N/S interface 
where the predicted feature in f(E) has its maximum 
amplitude. This location introduces a new difficulty in 
our experiment because the density of states (DOS) of 
the Ag wire near the N/S interface is strongly modified 
by proximity effect. Hence we must first determine the 
modified DOS at equilibrium before we measure f(E) 
under nonequilibrium conditions. 

The current- voltage characteristic of the probe tunnel 
junction is 



I(V ) 
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dEn M (E) J deP(e) 



P , , (10) 

en,T J J 

[f M (E)n As (E -eV-e)(l- f Ag (E - eV - e)) 
(1 - f Al (E))n As (E -eV + e)f Ag (E - eV + e)], 



where Rt is the normal state tunnel resistance, n Ag and 
n A \ are the normalized densities of states, and f Ag and 
/ai are the electron energy distribution functions on the 
Ag and Al sides of the tunnel junction, respectively. The 
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function P(e) characterizes the probability for an elec- 
tron to lose energy e to the resistive environment while 
tunneling across the oxide barrier, an effect known as 
"dynamical Coulomb blockade" M- P(e) was determined 
from equilibrium measurements at high magnetic field, 
where superconductivity is completely suppressed. De- 
tails of the fitting procedure used to extract P(e) were 
reported earlier* 1 ^ 
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FIG. 12: (color online) Expanded view of differential conduc- 
tance data for eV near A from the tunnel probe far away 
from superconducting reservoirs, for select magnetic fields. 
Symbols represent data while solid lines are best fits to BCS 
theory using a single value of the gap, A, and a depairing 
strength proportional to B 2 . Notice that the data at B = 
deviate significantly from the theory, whereas the data sets 
with B > 13 mT are fit well by the theory. 



FIG. 11: a) Differential conductance data and their best fit 
for the reference S/I/N tunnel junction at B — 13 mT and 
T = 40 niK. b) Blue line is the Al density of states, uai(E), 
used to produce the fit in part (a). (The black line shows the 
ideal BCS DOS without a magnetic field, for comparison.) 

c) The dl/dV data and their best fit in equilibrium for the 
tunnel probe on the sample shown in Fig.[JJ at B = 12.5 mT. 

d) Blue line is the riAg(E) used to produce the fit in part (c). 
The black line is a fit to the solution of the Usadel equation 
discussed in the text. 

Quantitative analysis of our tunneling data requires 
an accurate determination of the superconducting gap, 
A, hence we fabricated a second S/I/N tunnel junction 
simultaneously with the sample, but placed about 20/im 
away from it, and with the N side of the junction far 
from any superconductor. Tunneling spectroscopy mea- 
surements on this reference junction, shown in Fig. lllf a). 
were fit to Eq. (|10p with riAg independent of energy and 
with the standard BCS form for Uau to provide an accu- 
rate determination of A. 

Several of our tunnel junctions exhibited sharp anoma- 
lies in the conductance data for voltages close to the su- 
perconducting gap; however, these features disappeared 
with the application of a small magnetic field of B = 
12.5 mTj^ Figure [T^l shows dl/dV data for one partic- 
ular tunnel probe for different magnetic field strengths. 
Along with each data set are fits using the standard BCS 
form for uai(E) with a small depairing parameter pro- 
portional to B 2 , which has been shown to account well 
for applied magnetic fields^ Adding this term effectively 



rounds the DOS in the superconductor. Following the no- 
tation of Ref. |3^ we determined a depairing parameter of 
7 = T/A = 0.0020 for B = 12.5 mT and a supercon- 
ducting gap in the Al of A = 274 fj,eV. This rather large 
value for A was consistent across samples and is believed 
to be due to oxygen incorporated into thin, thermally- 
evaporated Al films. 40 ' 41 

With the form for uai(E) confirmed, it is possible to 
analyze the dl/dV data for the sample tunnel probe, 
which is in close proximity to superconducting reservoirs. 
Figure fTTT b) shows the dl/dV for this probe with an 
external field of 12.5 mT applied. The differences be- 
tween the dl/dV data from tunnel probes nearby or far 
away from superconducting reservoirs arise from changes 
in the DOS of the normal wire due to the proximity ef- 
fect. Rather than a flat DOS used to fit the data in 
Fig. [TlTa). the DOS in the Ag wire near the supercon- 
ducting reservoirs is modified as shown by the squares in 
Fig. fTlTd) . This shape was obtained by deconvolving the 
dl/dV data. 

Also shown in Fig. [TlTd) is the density of states of 
the Ag wire determined from a numerical calculation of 
the Usadel equation (solid line). This calculation re- 
quires knowledge of the sample dimensions, the gap in 
the superconducting reservoirs, and the Thouless energy. 
The sample dimensions were obtained from scanning elec- 
tron micrographs, such as the one shown in Fig. [TJ The 
gap in the superconducting reservoirs was found to be 
A w 150 /xeV. This value of A is much smaller than 
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the value in the Al tunnel probe because the reservoirs 
are much thicker than the tunnel probes (and presum- 
ably contain much less oxygen), and because they are 
close to a normal metal-superconductor bilayer. Finally, 
the Thouless energy was determined by fitting the criti- 
cal current vs. temperature data as discussed in Section 
IIV CI The value of Etk was then refined through self- 
consistent calculations involving both the finite probe 
size and the position dependent order parameter A in 
the superconducting reservoirs. 
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FIG. 13: a) Density of states of Ag wire at location of tunnel 
probe with different amounts of supercurrent flowing across 
the S/N/S junction. Solid and hollow circles represent tia 9 
for I a = 0.9/ c and I s = — 0.97 c , respectively, while solid line is 
for I s — 0. Inset: Theoretical results of injecting supercurrent 
into device. Solid line for I 3 — and dots for I s = ±.91 c 

When supercurrent flows through the S/N/S junction, 
the phase difference of the reservoirs, (f>, changes riAg(E). 
For this reason, dl/dV data were also taken with super- 
current flowing across the S/N/S junction. The resulting 
fits of riAg for I = 0.9I C and I = — 0.97 c , which are iden- 
tical to each other, are shown in Fig. [T3] along with one 
for I = 0. It is noteworthy that the change of shape 
(more narrow at low energies, broader at intermediate 
ones) is qualitatively consistent with theoretical calcula- 
tions shown in the inset. 

It was anticipated that applying a voltage to the nor- 
mal lead would not alter riAg(E) so that it would be 
possible to deconvolve the distribution function, f A g{E), 
for the system out of equilibrium. Figure [T?] shows that 
this assumption does not hold, as the best fits for applied 
voltages U = 22/iV and U = 63 \xV are poor. Only by 
using an altered UA g (E) was it possible to fit the data in 

Fig.rrn 

Changes in riAg(E) with increasing U are proba- 
bly due to a slight suppression of the gap in the 
superconducting electrodes, which are adjacent to 
superconductor/normal- metal bilayers. To estimate how 



FIG. 14: Differential conductance tunneling data taken with 
a voltage U applied between N to Si to drive the system out 
of equilibrium. Red lines are the best fits using the riAg data 
from Fig. [TJ] a) U = 22 fiV. b) U = 63 fj,V. The fits are 
unable to reproduce the data. 



riAg(E) changes, we used two different forms for the dis- 
tribution function fA g (E). First, we computed fA g (E) 
from Eqs. ©(a) and (b), which include proximity ef- 
fects due to the nearby superconducting reservoirs, but 
neglect inelastic scattering. Second, we solved the dif- 
fusive Boltzmann equation with collision integrals for 
electron-electron scattering, while neglecting supercon- 
ducting correlations. The two forms for the distribution 
functions are shown in Figs. I15l fa) and (b) for two dif- 
ferent values of U: 25 \iV and 63 fiV. Using those distri- 
bution functions, new densities of states were obtained 
by deconvolution of the dl/dV data and are shown as 
the symbols in Figs. fTBTc) and (d). Notice that the two 
different forms of the distribution function yield similar 
results for n^g. Figures [TBTc) and (d) also show %g ob- 
tained from equilibrium dl/dV data, as the solid lines. 
However, the resulting nA g does not obey the sum rule 
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that should be valid in all situations. We do not know 
what causes this discrepancy. 

Fortunately it is possible to extract information about 
f Ag (E) using a method that is relatively insensitive 
to the exact form of riAg(E), by taking advantage of 
a near-symmetry of the data with respect to the di- 
rection of Ig. The data shown in Fig. [T3] confirm 
our expectation that riAg(E,U — 0, Is) — riAg(E,U — 
0, — Is), a symmetry that also holds approximately for 
U 0. Hence, one can analyze the difference between 
two data sets with opposite directions of the supercur- 
rent, dI/dV(V, U, I s ) - dI/dV{V, U, -!<?), which will de- 
pend on the differences in the distribution functions, 
Sf Ag {E) ee f Ag {E,U,I s ) - f Ag (E,U,-I s ). The effect 
of analyzing the data with the wrong DOS for the Ag 
is greatly reduced in this case. The feature we seek in 
SfAg(E) is predicted to be odd in I s , hence it should be 
the only contribution to SfAg(E). Figure [TCT a) shows 
5f with U = 22 /iV and Is = 0.9I C , which exhibits the 
predicted feature that is antisymmetric in energy. The 
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FIG. 15: Calculated DOS using expected forms for distribu- 
tion functions, a) Distribution functions for U — 25 fiV. Blue 
line calculated by solving Eqs. © without collision integrals. 
Red line calculated by solving Boltzmann equation including 
collisions, but not including superconducting correlations, b) 
Distribution functions for U — 63 fiV. Below are the decon- 
volved forms for n Ag using the above distribution functions 
when c) U — 25 /iV and (d) U = 63 fj,V. In both, black lines 
represent U = (jV, for comparison. 



solid lines are the numerical solution to Eqs. ^j, with the 
parameters Eth and A obtained from the previous fits, 
with no additional fit parameters. The computed theory 
curves agree well with the experimental data. 

A further test of the robustness of the experimen- 
tal results is to compare the measured form of Sf Ag (E) 
when the signs of both U and 1$ are reversed, i.e. 
f Ag (E, -U, -I s ) - f Ag {E, -U, I s ). The results of this 
second measurement are shown superimposed on the first 
in Fig. HUa) . The agreement between the two data sets 
is excellent. 

Interestingly, applying a voltage U > 34 /iV brings this 
sample into its 7r-state. Figure UMb) shows 5f(E) data 
for U = 63 fiY and I s = 0.9J C . Compared to Fig. [TBJa), 
the sign of the low-energy feature in Sf(E) is reversed, 
demonstrating that the phase difference <fi, rather than 
the supercurrent, determines the sign of the new feature 
in f(E). 

The results of Fig. [TH] indicate that the supercurrent 
has a large effect on the electron energy distribution func- 
tion inside the normal metal. Such a mechanism has been 
utilized to explain^ the large thermopower measured in 
Andreev interferometers^ — systems with two normal- 
metal and two superconducting contacts. Our results 
confirm this mechanism and point out to new phenom- 
ena dependent on it, such as the large Peltier effect^: in 
linear response to the quasiparticle current (voltage), the 
supercurrent-induced change 5f translates into a change 
of the electron temperature and the sign of this change 




FIG. 16: (a) Sf(E) = f Ag (E, U, I s ) - f Ag (E, U, -I s ) for U = 
22 fiV and I s = 0.97 c . (b) Same quantity for U = 63 
and Is = 0.9/ c , where U > 34 /iV corresponds to the system 
being in the 7r-state. In both figures, a second data set (open 
circles) is shown with the signs of both Un and 7s reversed. 
Solid lines are numerical solutions to Eq. ©. 



(heating or cooling) depends on the relative sign of the 
supercurrent compared to the sign of the quasiparticle 
current. One can hence cool part of the structure by si- 
multaneously applying a quasiparticle current and a su- 
percurrent. 



VII. CONCLUSIONS 

Superconductor/normal metal hybrid systems exhibit 
a wealth of fascinating behaviors, starting with the prox- 
imity and Josephson effects. Driven out of equilib- 
rium, the possibilities increase, from the non-equilibrium 
7r-junction to the supercurrent-induced modification of 
f(E) discussed in the final section of this paper. All of 
these observations can be interpreted with two main con- 
cepts: the spectrum of the supercurrent Je{E) and the 
electron distribution function f(E). The latter can be 
tuned by applying voltages or changing the temperature 
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— the previous for example by applying a magnetic field. 
In Sees. IIII1 IIVI and IVII we showed in different schemes 
how the nonequilibrium f(E) changes the observed su- 
percurrent and how the supercurrent affects f(E). In 
Sec. |V] we discuss the modifications in Je{E) due to a 
magnetic field and the resulting changes in the supercur- 
rent. To our knowledge, the effect of a parallel magnetic 
field on the S/N/S critical current had not been explored 
in detail before. 

As discussed in Refs. Infill the Zeeman effect due to a 
magnetic field will cause analogous changes in the super- 
current as the nonequilibrium population of the super- 
current carrying states. This exact analogy is distorted 
on one hand due to the inelastic scattering changing the 
nonequilibrium distribution function, and on the other 
hand the orbital effect arising from the magnetic field. It 
remains an experimental challenge to show this analogy 
and combine the two effects in the case when the Zeeman 
effect dominates over the orbital effect. As we discuss in 
Sec. El the latter would require constructing extremely 
thin junctions. 

Recently there has been intense interest in the limit 
where a Josephson junction behaves as a coherent quan- 
tum system with one degree of freedom^ 3 - There is hope 
that Josephson junctions will someday provide the build- 
ing blocks for a quantum computer. In the meantime, 
we hope to have demonstrated that even in the classical 
regime, the Josephson junction is full of surprising new 
possibilities. 

While preparing this manuscript, we learned about re- 
cent related works^i where the magnetic field dependence 
of the S/N/S supercurrent was also studied. 
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APPENDIX A: USADEL EQUATION AND THE 
MAGNETIC FIELD 

The Usadel equation in a magnetic field can be writ- 
ten as follows, making use of the ^-parameterization 
G = cosh (9, F = e* x sinh# of the quasiclassical Green's 
functions) 4 *^^ 

v 2 

HD\7 2 = -2i(E + ah) sinh + (hT st + -2-)sinh20, 

2D (Af) 

V • (v s sinh 2 0) = , v s =D[V X - 2eA/h] . (A2) 

Here, D is the diffusion constant, A the vector poten- 
tial, r s f the spin-flip rate, v s the gauge-invariant su- 
perfluid velocity and h — -^gfi\B\ the Zeeman energy. 
The equation is to be solved separately for both spin 
configurations a = ±, assuming spin- independent mate- 
rial parameters. The spin-averaged spectral supercurrent 
jE — 5 Im[w s sinh 2 0\ a =+ + v s sinh 2 0\ a= J\/ D is obtained 



from the solutions and can be used to calculate the ob- 
servable supercurrent under various conditions. Below, 
we consider these equations in a wire that has an uniform 
cross-section S, and assume the boundary conditions 

X = ±0/2 , = 0o, at x = 0, L , (A3) 

h ■ v s = , h ■ V0 = , on dS , (A4) 

where OS is the boundary of S. These imply that we ne- 
glect details of the current distribution near the terminal - 
wire contact. 

When a magnetic field is applied to a wire, in addi- 
tion to the Zeeman splitting, the field induces circulating 
components to the supercurrent flowing in the wire (see 
Fig. fT7]) . These currents contribute to decoherence, for 
example reducing the magnitude of the critical current, 
and in the general case also prevent reducing Eqs. (|A1|) to 
one-dimensional equations in the direction parallel to the 
wire. For thin wires, however, the additional decoherence 
can be simply absorbed to the spin-flip parameter T s { in 
the one-dimensional Usadel equation and the vector po- 
tential can otherwise be neglected i 33 i 39 ' 44 i 45 Below, we 
show how this conclusion can be reached for an arbitrary 
orientation of the magnetic field, and that the results are 
consistent with the discussion in Section fVl 

Reducing Eq. (|A1|) to a one-dimensional equation is 
possible when the transverse dimensions d of the wire 
satisfy d <C L,l m , where L is the distance between the 
superconducting contacts and l m — yJh/eB a magnetic 
length scale. This is beca use varies on the length scales 
of l m and I e = y/hD/E ~ L when considering energies 
E ~ hD/L 2 relevant for the supercurrent. For perpen- 
dicular fields it is also possible to directly choose a proper 
London gauge where x varies slowly in the transverse di- 
rection and vs oc A. 

Since l rn ~ 80 nm for B ~ 0.1 T, in the experimen- 
tally interesting situation we have w < l m <C Ie ~ L. 
To handle the details of the problem in this case, we 
apply perturbation theory in the parameter A = d/L. 
We choose a coordinate system such that x is the coor- 
dinate parallel to the wire and y and z correspond to 
transverse directions, and fix a convenient gauge A = 
(B y z — B z y, —B x z, 0) in which the vector potential is in- 
dependent of x. Finally, we rewrite Eqs. (|A1[) in the 
dimensionless variables x = xj L, y — y/XL, z — z/XL, 
B = eL 2 H/hX and substitute in the regular series expan- 
sion = O + X0i + X 2 2 + . . ., x = Xo + Xxi + X 2 X 2 + ■■■■ 

Requiring the equations corresponding to orders A -2 , 
A -1 and A of expansion to be separately satisfied, we 
first find that the variables O , 0i and xo are independent 
of y and z. We also find that the first-order response 
f>X = ^Xi is given by 

Vi^x-O, h- V±5x\dS = n-2eA/h. (A5) 

Here, dS is the boundary of the cross-section of the wire, 
h its outward normal vector, and the operator Vi con- 
sists of the transverse components of the gradient. This 
x-independent result applies in the central parts of the 
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wire, away from boundary layers near the ends of the 
wire. Finally, after averaging the equations of order A 
across the cross section S of the wire, we arrive at the 
result 

d% = -2i(E + ah) sinh + (hT si + 7) sinh 29 , (A6) 



^(sinh 2 6 d x xo) = 0, 



(A7) 



HD 



= 2S I d V dz O^XO + V±(5x - 2eA/hy , (A8) 
J s 

This shows how the effective spin-flip parameter is mod- 
ified by the applied field B. 

The additional decoherence (|A8[) depends on the di- 
rection of the field and the cross section of the wire. For 
a wire with a circular cross section of radius R, we note 
that Eq. (|A5[) has the exact solution 5\ = —eB x yz/h. 
This results to 

7 = \hD{d xXQ ? + ?R (±R 2 B 2 X + B?Bl + R 2 B 2 Z 

(A9) 

For wires with a rectangular cross section, we cannot 
solve Eq. (|A5|) analytically. However, a variational solu- 
tion is still possible: we can expand 8\ in polynomials 
of y, z to orders n < 3 and project Eq. (|A5|) onto this 
function basis. From this procedure, we find 



2eB r 



h dy + d 2 



yz, 



7 - -hD(d xX o) 2 + {w 2 yz B 2 x + d 2 z B 2 y + d 2 y B 2 z ) 



d 2 d 2 

W ~ 

yz dl + dl 



(A10) 

(All) 
(A12) 



Approximations using higher-order basis produce only 
slight improvements in accuracy to 7. 

We now note that the contribution of the magnetic 
field to the decoherence rate 7 is of the form e 2 Dd 2 B 2 /h 
for all directions of the field, where d is proportional to 
some transverse dimension of the wire. Comparing this 
to the energy scale = hD/L 2 of the one-dimensional 
Usadel equation (|A6|) , we find a dimcnsionless parameter 
(eBLd/h) 2 oc (<&/<I>o) 2 that determines how much the 
magnetic field suppresses coherence. Here, the flux $ 
corresponds to an area Lx d, which is in agreement with 
the discussion in Section fVl 

Finally, note that above we neglected the screen- 
ing of the magnetic field by the induced supercurrents. 




FIG. 17: Supercurrent flow induced by a magnetic field B = 
(B x , By, B z ) oc (3, 1, 2) in a thin rectangular wire. The arrows 
indicate the magnitude and direction of the superfluid velocity 
vs- Fourth-order variational solution for \ i s used here, see 
text. 



where d y and d z are the width and thickness of the wire. 
For orders n < 4 we obtain instead 



d 2 d 2 

~2 _ y z 



266d 2 d 2 



1054 + 1500^ + 105d 4 



(A13) 



However, this should not be important in the experi- 
mental case, as the Josephson screening length Aj = 
yjhd 2 /2efioI c L > 200 nm is larger than the width of 
the junction, and the aluminum terminals are sufficiently 
thin as to produce only small screening. 
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